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Abstract. A reference potential approach to the one-dimensional quantum- 
mechanical inverse problem is developed. All spectral characteristics of the 
system, including its discrete energy spectrum, the full energy dependence of 
the phase shift, and the Jost function, are expected to be known. The techni- 
cally most complicated task in ascertaining the potential, solution of a relevant 
integral equation, has been decomposed into two relatively independent prob- 
lems. First, one uses Krein method to calculate an auxiliary potential with 
exactly the same spectral density as the initial reference potential, but with 
no bound states. Thereafter, using Gelfand-Levitan method, it is possible to 
introduce, one by one, all bound states, along with calculating another auxil- 
iary potential of the same spectral density at each step. For the system under 
study (diatomic xenon molecule), the kernel of the Krein integral equation can 
be accurately ascertained with the help of solely analytic means. At small 
distances the calculated auxiliary potential with no bound states practically 
coincides with the initial reference potential, which is in full agreement with 
general theoretical considerations. Several possibilities of solving the Krein 
equation are proposed and the prospects of further research discussed. 



1. Introduction 

In this paper, a recently proposed reference potential approach to the one- 
dimensional quantum mechanical inverse problem [P is further developed. Let 
us briefly recall that the starting idea is to choose a suitable reference potential 
for the system. For this fixed potential, it is always possible to calculate all its 
spectral characteristics, which means that in this artificial way one can obtain the 
complete set of information (otherwise inaccessible) needed to uniquely solve the 
inverse problem: 1) full energy spectrum of the bound states; 2) full energy de- 
pendence of the phase shift for the scattering states; 3) norming constants of the 
regular energy eigenfunctions for all bound states. Of course, there is no need to 
regain the potential which is already known by definition. However, the quantities 
related to the reference potential can be used as zeroth approximations to the real 
spectral characteristics of the system. A further step might be, for example, to 
calculate another potential whose discrete energy levels would exactly fit with their 
actually observed values, so that at least in this sense the new potential would be 
more realistic than the initial reference potential. 

To illustrate the method, as previously, Xe2 molecule is under study, and the 
same three-component exactly solvable reference potential is used. We have already 
ascertained the full energy dependence of the phase shift for the scattering states 
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PQ and demonstrated its excellent agreement with the celebrated Levinson theorem 
In addition, we are provided with full knowledge of the Jost function, which 
is the most important spectral characteristic of the system. Thus, we are prepared 
to attack the most serious computational-technical problem, solution of an integral 
equation, which would enable to uniquely ascertain the potential. For this purpose 
a combined approach is used. In Section 2 we recall some useful properties of the 
Gelfand-Levitan method [2], which enable to separate the main problem into two 
independent parts: 1) calculating an auxiliary potential Vo(r) with exactly the same 
spectral density for positive energies as the reference potential V(r), but with no 
bound states; 2) calculating a series of auxiliary potentials Vk(r), k = 1,2, ...,n 
(i.e., adding, one by one, all bound states, and keeping their norming constants), 
also having the same spectral density for positive energies. The second step is, in 
fact, less complicated and is only briefly discussed in this paper. Our main goal is 
to ascertain the auxiliary potential Vo(r) with no bound states, starting from the 
known Jost function for positive energies. This is the subject of Section 3 where 
the Krein method 0| E] will be used. Some concluding remarks and a discussion of 
the further perspectives of the method form the content of Section 4. 



2. Step-by-step use of Gelfand-Levitan method 

First, let us recall general solution scheme in the frame of Gelfand-Levitan 
method |Hj (see, e.g., [Hj for more details). One's aim is to solve the integral equation 



(1) 

whose kernel 
(2) 



K{r,r') + G(r,r') + / K(r, s)G(s,r')ds = 0, 



G(r,r') = 



sin (kr) • sin (kr') 
k~ 2 



da, 



is determined by the quantity da = 
the difference of two spectral densities 



dp 2 (E) d Pl {E) 



dE 



dE 



dE, which contains the 



(3) 



dp{E) _ J n- 1 ^\F(E)\ 2 ,E>0, 
~dE~ ~ 1 £ C n S{E - E n ), E<0, 



one of them ( - \ g - ) being related to a known potential (Vi(r)). In this sense 
dE 

Gelfand-Levitan method also represents a reference potential approach. The sim- 
plest possibility is to take V\(r) = 0, and correspondingly, ^ = . Then 

dE 7T 

00 

(4) G(r, r') — — / sin (kr) ■ sin (kr 1 ) g(k)dk + > — % sinh Oy n r) sinh Oy„r') , 
n J ' „ 4 7 2 
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where 7^ 



2mE„ 



E n being the bound levels and C n , their related norming 



constants. The characteristic function 
(5) g(k) ee 



\F(k)Y 



is determined by the modulus of the Jost function F(k). Thus, in principle, one 
can solve Eq. (1) and then calculate the potential 



(6) 



V(r) = 2C-^K(r,r) 1 C=^-. 

dr 2m 



This scheme might seem simple, but its actual realization stumbles upon serious 
computational-technical difficulties. Now, let us assume that we somehow managed 
to ascertain an auxiliary potential Vo(r) with exactly the same spectral density for 
positive energies as the desired potential V(r) 7 but having no bound states. In this 
case one easily finds that |B] 



(7) 



G(r, r') = C i<Po {ijj,r) <^ {il 3 ,r') 



where (po (i"fj, r) are the regular solutions (NB! not the real eigenfunctions!) related 
to the auxiliary potential Vo(r). Next, let us introduce another auxiliary potential 
V\{r) with just one bound eigenvalue E\ (note that numeration of the levels here 
starts from the highest-energy one, and therefore, E n corresponds to the zeroth 
level), and again, with exactly the same spectral density for positive energies (|F(fc)| 
remains the same). Thereafter, one introduces an auxiliary potential Viir) with two 
levels {Ex and E 2 ), etc., until he comes to the desired potential V n {r) = V(r). 
It can be proved that 

(8) Vj(r) - V (r) = -2C {ln[det Cjir)]}" , j = l,2,...,n 

and the corresponding regular solution 



(9) 



tfj (k, r) = det 



Ci(r) *i(r) 
(3j(k,r) (fio(k,r) 



[detQ(r)]- 1 



Here, a j x j matrix Cj(r) = 1 + f R^'(s)ds (I is the unit matrix) with the elements 



of R,W) being 

(10) R { i2, = Cm (ill,r) ifo (ilfm,r) , (Z,m = 1,2, 

a column vector (do not confuse the norming constants Ci with the matrices Cj{r)) 



(11) 



Mr) 



( \ 

V Cjipo {il 3 ,r) J 
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and a row vector 
(12) 

(r r r 

J <A) («7ii s ) ¥>o (fc, s) ds y (fi (hz, s) <fo {k, s) ds ... J tp (ijj,s) <p (k, s) ds 


have been introduced. 

To be more specific, let us examine the simplest case of just one bound state 
(j = 1). Then Cj(r) reduces to scalar and one gets 



(13) 



V 1 (r) = V (r)-2C{ln 



1 + Ci / tpKiji, s)ds 



As can be shown the regular solution ip% (171, r) related to V\{r) (i.e., the real 
confined eigenfunction) reads 



(14) 



fi (*7i>r) 



¥0 (iji, r ) 



1 + Ci / (fl {i-f!,s)ds 




Let us prove that Ci, as needed, is the norming constant of ipi (iji, r). Indeed, 
slightly rearranging Eq. (14): 1 + C\ f b« (171, s)ds = ; differentiat- 

ing both sides: C l( ^ r) = jgoyi_jgojgi ^ i £?i (171 , r) = - (^) ~ 

— tp[ = — ( — j , and integrating both sides of the latter equation, one gets 

</?o V^o/ 

00 

j C\(fi («7i,r) <ir = 1, because any regular solution ip(r) w r as r — > 0, and thcrc- 


fore, lim r ^o ( — ) = 1- 
\<PoJ 

Now, let us analyze a more complicated case of two bound states (j = 1). Ac- 
cording to Eqs. (8) and (13), 
(15) 



V 2 (r)-V 1 (r) = V 2 {r)-V Q (r)-[V 1 (r) - V (r)} = -2C { In 

'Cu(r) C 12 (r) 



detC 2 (r) 



l + Cifipl (iji,a)ds 




Here, the 2x2 matrix C 2 {r) 



C 2 i(r) C 22 (r) 



has the elements 



(16) C n (r) = l + Ci J (p (i'y-L, s) ds, C 12 (r) = Ci J ip Q (171, s) (p {i^ 2 , s) ds 


r 

C 21 (r) = ^C ia (r), C 22 (r) =1 + C 2 f ^{i l2 ,s)ds, 
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and therefore, the argument of the logarithm in Eq. (15) reads 

n 2 



C\C 2 



(17) 



fir) = C 22 - 



f(p («7i,s)v? {i"f2,s)ds 
.0 



1 + Cijifl (iji,s)ds 



1 + C 2 < / (fil {i-j2,s)ds 



Ci[(3i2(r)Y 



1 + (ij 1: s)ds 



where /?i 2 (r) = J" <^ (i7i> s) ^0 («72, s) oLs. 


As can be easily proved, 



f(r) = l + C 2 / <pl{i l2 ,s)ds 



(18) 



Indeed, using Eq. (14) and a more general formula [H] 



/ <^ («7i7s)^o (k,s) ds 



<pi (k, r) = ip (k, r) - dipo (iji,r) x — 



l + Ciftfl (iji,s)ds 




and taking k = 272, one gets 

<Pi [il2,r) = (f [ij 2 ,r) - Cnp x (i<ji , r) f} 12 (r) , 



i.e., 



<fl [il2,r) = (pi (i72,r) + [Ci^i (iji,r) 0is{r)] - 2C\ip Q (ij 2 ,r)<pi (iji,r) 0u{r) 



tpo («7i,r) 

which, after integrating from to r, and taking account of Eq. (14), proves the 
claim. 

Thus, in full analogy with Eq. (13), one can write 



(19) V 2 (r) = V x {r) -2C { In 

and more generally, 



r 

1 + C 2 J <pl (ij 2 ,s)ds 



(20) V h (r) = V k - 1 {r)-2C{hx 



l + C k / tp k _ 1 (ij k ,s)ds 



k = 1,2, 71. 



This way one gets rid of inconvenient matrix equations and can move towards the 
goal step-by-step, introducing one bound state at a time. 

The described procedure can be reverted 0|H|, i-e., starting from the potential 
Vfe(r), one can construct the auxiliary potential Vk-i[r), removing the zeroth level 
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(Ek according to the numeration used). The new potential reads 
(21) 



V k _ 1 (r) = V k (r)+2C 



where "J/^ (r) is the eigenfunction of the zeroth level, not necessarily normalized 
(note that the norming constant is absent here). One can remove, one by one, 
all bound states until he comes to the potential Vo(r) with no bound states. An 
illustration to Eq. (21) can be seen in Fig. 1. 

Thus, we have demonstrated that there is no need to solve Eq. (1) all at once. It 
is probably much easier to first solve this equation for the auxiliary potential Vo (r) 
and then, step-by-step, introduce the bound states as described above. The Jost 
functions F n (k) and Fo(k) related to the potentials V n (r) and Vo(r), respectively, 
are connected by a simple formula 



(*«(r) 



ds 



(k = n, Ti—l, 



(22) 



F„(k) = F (k) Yl 



k- 



*7i 



Hi 



which demonstrates that these potentials, indeed, have exactly the same spectral 
densities for positive energies . On the other hand, Eq. (22) illustrates the general 
rule that any zero of the Jost function corresponds to a bound state. 

To end this section, let us recall important asymptotic formulas (see [§] for 
details) 



(23) 



V n (r)-V (r) 



-AC 



r, r 



0. 



(24) 



2C 

V„(r) - V {r) « -— (2 7 i) 5 exp(-2 7l r) I r 



From Eq. (23) one can infer that in the immediate vicinity of the zero point r = 
the potentials V n (r) and Vo(r) practically coincide. Indeed, all norming constants 



Ci = 



f<Pj (hj,r) dr 



are extremely small quantities, since they are related to 

oo, but 



the regular solutions proportional to r as r — > 0, and vanishing as r 
achieving very large absolute values between these asymptotic regions. 



3. Krein method 

Thus, in view of the results of previous section, our main goal is to accurately 
ascertain the auxiliary potential Vo(r) with no bound states, provided that its 

spectral density ^ — - = -k^ x \/E \F[E)\ 2 is known. This is indeed the case, 
dE 

since we have carefully calculated the modulus of the Jost function for the reference 
potential V(r) £Q, and the whole idea is to use the same quantity (\F(E)\) for the 
auxiliary potential Vo(r) as well. Since there are no bound states, the kernel of the 
Gelfand-Levitan equation, according to Eq. (4), becomes analytically very simple: 

(25) G(r,r') = H(r-r') - H{r + /), 
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where a new function 



(26) H(r) = tt- 1 J g{k) cos(kr)dk, 

o 

with g{k) given by Eq. (5), has been introduced. The function H{r) = H(—r) is 
very important for the further treatment, wherefore a special designation, Krein 
-ff-function, will be used for this item henceforward. A motivation for this name 
stems from a very useful method for solving the inverse problem, which has been 
elaborated by Krein [5]. This method is based on a Fredholm-type integral 
equation 

(27) r 2r (r') + H(r') + J T 2r (s)H{s - r')ds = 0. 

o 

It can be shown that the solutions of Gelfand-Levitan and Krein equations are 
linked by a simple formula (cf. with Eq. (25)) 

(28) K(r, r') = T 2r (r - r') - T 2r (r + r'), 
and the desired potential reads 

(29) Vo(r)=A c([G(x)} 2 -^±\, x^2r, 



dx 

where G(x) = T2 r (2r). A special notation, Krein G- function, will used for this 
quantity henceforward. 

Without any doubt, Krein method is well suited for our purposes, but there is 
still a lot of analytical and computational-technical work to do. First, we have to 
accurately ascertain the Krein H- function, which, according to Eq. (26), is simply 
the Fourier cosine transform of the characteristic function g{k) whose full energy 
dependence has already been ascertained [I] . One might think that the H- function 
can be determined using the well-known fast Fourier transform technique. This, 
however, is an erroneous view, because (as we demonstrate below) the 77-function 
has to be calculated in a wide distance range with very small step, to ensure the 
correct asymptotic behavior of the resulting potential. Fortunately, the problem 
can be solved accurately and quite easily with the help of solely analytic means. 

3.1. Calculation of iJ-function. The function g(k) is shown in Fig. 2. In addi- 
tion to the overall curve, some characteristic slices can be seen in the insets. There 
exists a range k € (0, k\) where g(k) = — 1 (with high accuracy), and for this range 
one immediately gets the relevant component of the ii-function 

sin(fe x r) 

H o( r ) = • 

ixr 

Now, let us see how the g(k) curve passes through the "critical" region around 
ka = y/V(0)/C. As it happens (see the upper inset in Fig. 2), in a narrow range 
k £ (ki, k 2 ) (fco is also located within this range) g(k) + l can be nicely approximated 
by a sum of several Gaussians: 



(30) = -! + $>. 



j exp 
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where the parameters Oj, bj and kj can be determined from a least-squares fit. For 
the reference potential examined here, four such components have been introduced, 
and the value fc 2 = 19230 A -1 has been chosen. Thus, in the above formula one 
can replace k\ with fc 2 (due to -1 in Eq. (30)), i.e., 

sin(/c2'™) 



(31) 



fli(r) 



while the sum of Gaussians gives another component of the Krein H-function, 
H 2(r) = J2 H 2\ r )> where all constituents H { 2 3] (r) can be ascertained analytically. 

3 

k — k- 

Indeed, one can introduce a new independent variable x = — _ ■ and calculate 
(32) 



H V\ r ) = y exp(-a; 2 ) cos (V^bjxr + Ayr) da; 



a A exp(-i& 2 r 2 )x 



where X; 



|cos(fcjr) [Reerf(y 2 ) - Reerf(t/i)] - sin(fc J r) [Imcrf(y 2 ) - Imerf(yi)] j , 
ki - k 



V2bj 



and yi 



xi 



ibjr 



— — 1,2 and i is the imaginary unit). 
V2 



The error function, for any complex argument z, can be evaluated in terms of the 
confluent hypergeometric functions 0| 



(33) 

where ^(a, c; x) = 1 ■ 
is more convenient: 



Brf(*) = ^Z* Q,| 



ax a(a + l)x 2 



(34) 



erf(z) = 1 



l!c 2!c(c + l) 

exp(— z 2 ) 



1 



1 

2^ 



For large arguments another expression 
1-3 1-3-5 



(2z 2 ) 2 (2z 2 ) 3 



Consequently, any constituent of the Krein TJ-function expressed by Eq. (32) can 
be easily ascertained with any desired accuracy. 

Next one can introduce an arbitrary (but still reasonable) boundary point k$, 
and approximate g(k) in the range k £ (fc 2 , £3) as follows: 



(35) 



.9 



(fc) = a,j cxp —bj(k — kj) 



This brings along another component of the iJ-function H^{r) — ^H^\r) with 



constituents 



(36) 



Hg'(r) = -*- / exp -bjik-kj) cos (At) dk 



ajexpibjkj, 



,{b) 



x {exp(— bjk$) [r sin (k^r) — bj cos (k^r)] — exp(— bjkz) [r sin (fc 2 r) — bj cos (fc 2 r)]} . 

Introducing new suitable boundary points etc., the approximation of g(k) 

in the form of Eq. (35) can be continued until the conventional starting point k a of 
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the asymptotic region (see the treatment below). The number of these boundary 
points, as well as the number of exponents in any particular interval, is, of course, 
a subject for probes and trials. We have introduced four such intervals and used 
a different three-exponent approximation in any of them. As has been carefully 
checked, this ensures the accuracy of at least 6 significant digits for the calculated 
Krein iJ-function in the whole physical domain. The parameters of all components 
are given in Table 1. 

For the remaining part of the fc-space the asymptotic formula 

07) = % + % + k> ka , 

can be used, where P 

V(0) 2[V(0)] 2 -CV"(0) 

(38) fl2 ___ a4 _ __ . 

The coefficient a^, as well as the coefficients for higher-order terms, can also be 
calculated in terms of the reference potential and its derivatives, but the resulting 
expressions are rather complicated and inconvenient for practical use. Instead, we 

Q,Q 

only introduced just one additional term, — ^, and determined the coefficient ae 

k° 

from the general demand (see Eq. (23)) that the potentials Vo(r) and V(r) should 
coincide as r — > 0. Such physically well motivated choice of is indeed possible, 
as will be explained below. Thus, within this approximation, the asymptotic part 
of g{k) reads 

(39) g{k) + + k > kai 
where 

(40) bi = -2a 2 , b 2 = -2(a 4 - a§), b 3 = -2(a 6 - 2a 2 a4 + ~a%). 



3 



The relevant asymptotic component of the 77-function becomes 

oo 



sin t 

Since the expression for H a (r) will contain the sine integral Si(x a ) = J dt, 

o t 

let us recall a useful formula 9 

(41) Si(a;) = \~\ exp(-ia?)*(l, 1; ix) + % - exp(za;)*(l, 1; -ix), 

where i is the imaginary unit, and the function ^(a, c; z) is a particular solution of 
the confluent hypergeometric equation introduced by Tricomi. For a large argument 
it can be evaluated from the asymptotic series 

JV 

(42) *(a,c;z) = z~ 

n=0 



-ay-' (fl)n(a- C + 1) T , 

^ n\{-z) n 



where (a)„ = T(a + n)/T(a) = a(a + l)(a + 2). ..(a + n — 1) is the Pochhammer 
symbol, and N must not be too large. If Eq. (41) is usable, i.e., in the case of 
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sufficiently large r, one gets the following expression: 

/ aq\ it i s fbiAx b 2 A 2 b 3 A 3 \ sin(x a ) 

( 43 ) *H a (r)= [——-—— + 



where 



l\k a 3!fc3 h\kl J x 
b x Bi b 2 B 2 b 3 B 3 \ cos(x a ) 



llk a 3!fc3 5!fc5 



, x a — k a r >> 1, 



j=0 \ Xa ) j=0 \ Xa ) 

If Eq. (41) cannot be used, one can apply the universal expansion 

(44) Si(x „ ) = x „_M! + ££_..., 

to get another formula 



b\r b 2 r 3 b 3 r 5 



(45) ^(r) = ( ^_ _ _i_ + _i_ ) &(*„)-- 



7T 



2 



+ 



6xXo 6 2 Xx , b 3 X 2 \ fb 2 Y bsYi. . 

cos(x a ) - r — ? - — T sm(x a ), 



(46) X, = £(-1)^ [2(i - j)]! (x a )« , = ^(-1)^' [2(i - i) + 1]! (z a ) 2j . 

From Eq. (45) one can infer that near the zero point r = 

H (r) ra tt" 1 f &1 I &2 I — ^ - — 
\fc a 3fc3 5fc a / 2 

and consequently according to Eqs. (38) and (40), the derivative at zero point 

[i? a (0)]' = — tJ- = ~^7T • Let us prove that the same relation holds for the overall 
z 40 

Krein iJ-function. Indeed, as r — > 0, one can always choose a value for fc a which is 
large enough, so that Eq. (39) can be used, but on the other hand, small enough, 
so that sin(fcr) w kr, if k € (0, k a \. Integrating by parts (note that 17(00) = 0), one 
gets from Eq. (26) 

k a 00 

nH(r -» 0) = -Jg'{k)kdk + 1 J + sin (fcr)&. 

fc a 

Only the second term gives contribution to H'(Q), and after few elementary trans- 
formations one comes to the desired result 

00 

(47) H (r^0) = - b -^+TT- 1 J g(k)dk, 



which proves that 

<48) - -h = m. 

Thus, in the case examined here, the Krein //-function can be, indeed, ascer- 
tained analytically. To get the overall //-function, one just sums the components 
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described above: H(r) = Hj( r )- The result for different distance regions is shown 

3 

in Figs. 3 and 4. As can be seen, H (r) is a rapidly oscillating function with decay- 
ing amplitude. The period of oscillations stabilizes quite rapidly and remains very 

close to the characteristic value L = — (see Eq. (31)), while the amplitude slowly 

approaches zero as r — > oo. 

3.2. Solution of Krein equation. Before setting about the main task, we have 
to fix the coefficient 63 in Eq. (39). To this end, let us take into consideration 
that according to Eq. (27), G(0) = —H(Q). Also, as repeatedly mentioned (see 
Eq. (23)), the potentials Vo(r) and V(r) should coincide as r — > 0. One therefore 
can replace Vo(r) with V(r) in Eq. (29), when studying the region very close to 
r = 0. In addition, the Krein G-function there can be approximated by a quadratic 
function: G(x) = a + bx + ex 2 (x — 2r), and consequently, G'(x) — b + 2cx. 
In this region a pseudo-Morse approximation for the reference potential is used: 
V(r) — Vo + Aq exp(— aox) — \J AqEq exp(— cxqx/2). Here, Aq = Dq exp(2aoro) and 
eo = Dq/4 (see pQ for more details). Therefore, the parameters a, b, c can be 
ascertained directly from Eq. (29). Indeed, one uses the relations N% = a 2 — b = 

^N 2 ^c-ab= a ^-^ /2 \N^b 2 + 2ac= al{AA ^^\ 
AC 8C ' 32C 

then solves the equation (a 2 — N-A + 2a [a(a 2 — N\) + N 2 ] — N 3 = to ascertain 

the parameter a = G(0) = —H(Q), and thereafter finds b and c. Having fixed 

these parameters, and also the value of k a , the coefficient 63 can be quite easily 

determined on a trial-by-trial basis. This way the value 63 = -5.883044-10 24 (A -6 ) 

has been fixed, which corresponds to k a = 75000 A -1 . 

Now, let us proceed with solution of Eq. (27). This equation can be discreticized 

using, for example, a four-point quadrature rule |10j 

(k+3)h 

(49) J f( x )dx=^f(kh) + ^f([k + l]h) + ^f([k + 2}h) + ^f([k + 3}h), 

kh 

which is exact for f(x) a cubic polynomial. Applying this formula to Eq. (27), one 
gets the following system of linear equations: 

(50) (1 + AH ) r 3n , + 3 (Ai^IVi + AH 2 r 3n , 2 ) + 2AH 3 T 3n>3 + 

+ 3 {AH^Y^n^ + Ai/ 5 r 3nj5 ) + 2AH e T 3n ^ + ... + AH 3n T 3n ^ n = — Hq 

Ai/_ 1 r 3n , +(i + 3AH ) r3„ ) i+3Atf 1 r3„ ) 2+2Atf 2 r 3 „,3+...+A J ff3n-ir3„ ) 3„ - -H t 

Afr_ 2 r 3Bj o+3Afr_ 1 r 3llil +(i + 3AH ) r 3n , 2 +... = -h 2 
A J ff_ 3 r 3 „,o+3(A J H'_ 2 r3 n ,i+A J ff_ 1 r3 n> 2)+(i + 2AH ) r 3n , 3 +... = -h 3 

AH- 3n T 3nt o + 3(A77i_3 n r 3ni i + A7J 2 _3nr3„^2) + (1 + AHq) r3 n> 3„ + ... = —H 3n 
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or, in a more compact form, 

(51) ^3n,k + A • (H-kTsnfi + 3Hi^k^3n,l + 5H 2 ~k^3n,2 + 2H 3 ^k^3n,3 + ••• 

+ H 3n ^k^3n,3n) = ~Hk, k = 0, 1, 2, 3n. 

Here A = — , H k = H_ k = H(kh), and r 3 „ ifc = T 2r (kh) (see Eq. (27)). 

o 

For any argument x = 2r = 3nh, Eq. (50) can be solved with the help of Gauss- 
ian elimination procedure, which is appropriate here, since only the last element 
of the solution vector, r3„.3„ = r2, (2r) = G(x) is actually needed to calculate 
the potential according to Eq. (29). Alternatively, one may rewrite Eq. (51) in a 
matrix form 

(52) (I + A-U)-G = -H, 
where I denotes (3n + 1) x (3n + 1) unit matrix, 

/ \ / 



G = 



^3n,0 
^3n,l 

V r 3 „ 3n / 



H = 



Hi 
V H 3n J 



and 



U = 



Hq 3Hi 3H 2 2i/ 3 ... 3Hs n -i H 3n 
Hi 3Hq .3Hi 2H 2 •••3i?3 7 j_2 Hs n -i 



\ H 3n 3H 3n _i 3H 3n _ 2 2_ff 3 „_ 3 ... 3Hi H ) 
The solution of Eq. (52) reads 

(53) G = - {I + A • Uy 1 -H = -H + A- U- H-A 2 -U 2 -H + 

and therefore, 

(54) 

^3n,3n = —H 3n + A ■ (H^nHo + 3H 3n -lH\ + 3H 3n - 2 H 2 + 2i?3„_3_ff3 + ... + HoH 3n ) 

— A 2 ■ g i gjH?, rl - l H\ l _ : j\Hj + ... 

The coefficients gi (here and henceforward) are defined as follows: gi = 1, if i = 
or i = m = 3n; gi = 2, if i = 3k and k = 1, 2, ...,n — 1; gi — 3 in any other case. 

According to Eqs. (6) and (28), V(r) = Acl^A dT *^ 

I dx 

this with Eq. (29), one gets an important relation 



dx 



. Comparing 



(55) 



£M = [G{x)]\ (x = 2r) 
dx 



Thus, instead of searching for the last element of the solution vector G, one may 
calculate its first element T m0 (m = 3n), which in some sense is more convenient. 
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Indeed, using such an approach, we can make use of the results of previous calcu- 
lations. Namely, as can be proved 

(56) r m+3 , = T m0 + A • [H m + 3H^ +1 + 3H^ +2 + H m+S ) - A 2 • S m + ... 

where 
(57) 

771—1 

S m = 2 gjHj(H m H m —j + 3H m +iH m +i-j + 3fl m -|_2-Hm+2— j + H m+3 H m+3 ^j)+ 
j=o 

Ho(3H^ + 9H 2 l+1 + 9H? n+2 + H 2 l+3 ) + 6Hi(H m+3 H m+2 + 3H m+2 H m+1 + 2H m+ iH, 

+6H 2 (H m+3 H m+ i + 2H m+2 H m ) + 4H 3 H m+3 H m , m = 0, 1, 2, ... 
From Eqs. (55) and (56) one gets another formula 

(58) 8 [G(x)} 2 =Hl + 3ff 2 i+1 + W 2 m+2 + H? n+3 - A • S m + x = mh, 

which may prove very useful, if the higher order terms can be ignored. 

In Fig. 5 one can see the calculated Krein G-function in the range from to 
1CP 6 A. Eq. (50) has been solved by Gaussian elimination, and an extremely 
small step h = 1CP 9 was used to ensure high accuracy. For comparison, another 
curve (actually, almost straight line) is depicted, which exactly corresponds to the 
reference potential V(r), i.e., it is the solution of Eq. (29) interpreted as Riccati 
equation: 

(59) ^M = [ G (,)] 2 -^, G (0) = ^(0). 

The curves practically coincide, which is a clear evidence of the validity of the 
approach. Indeed, this seemingly trivial calculation is, in fact, very sensitive to 
even minor inaccuracies in calculating the Krein 77-function, which would result in 
undesired and unphysical discrepancies, e.g., oscillations, of the G-function. Since 
no such discrepancies are seen at small distances, one may expect that both the 
Jost function and the TJ-function have been ascertained quite correctly. This in 
turn gives ground to hope that the potential can be correctly ascertained at longer 
distances as well. This, however, is not at all an easy task, and was not attempted 
here. 



4. Conclusion 

In this paper, as well as in the previous one PP, the treatment was more concen- 
trated on principles rather than the methods and techniques of computation. Wc 
demonstrated that the proposed reference potential approach enables one to ac- 
curately ascertain the important spectral characteristics, needed to uniquely solve 
the quantum-mechanical inverse problem. This way, one may get reasonable initial 
guesses to the real spectral characteristics of the system, which can be used, for 
example, to calculate another potential (Bargmann potential) whose Jost function 
differs from the initial one only by a rational factor. Let us briefly analyze the 
simplest case when this factor reads 

k — ia k — ia k + ia 

where both, a and b, are real and positive. In fact, it means that a discrete eigen- 
value — b 2 is replaced by —a 2 . As we see, Eq. (60) involves two operations. One of 
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them, introducing a new bound state E — —a 2 , can be performed as described in 

k + ia 

Section 2, while the second factor, -, causes an additional deformation of the 

k + ib 

initial potential [H] 

(61) AV(r) = V 2 (r) - V l{ r) = -2C {in W M« r ±^l r l }" , 

where the symbol W denotes Wronskian determinant and fi(ia, r) is the Jost solu- 
tion of the Schrodinger equation (/i(ia, r) — > cxp(— ar) as r — > oo). Both fi(ia,r) 
and the regular solution (pi(ib,r) are related to the potential Vj_(r). Thus, if —a 2 
is expected to be more realistic eigenvalue than — b 2 , one may hope that V^(r) is 
more realistic potential than V±(r) as well. 

The proposed approach is based on various analytic procedures. For example, we 
used piecewise analytic approximation for the characteristic function g(k) (see Eq. 
(5) and Fig. 2), and derived simple formulas for relevant constituents of the Krein 
H- function. In addition, we proved some general relations, Eqs. (47) and (48), 
regarding iJ-function near the zero point r = 0, and fixed an appropriate value 

oo 

for G(0) = — -ff(O) = 7T _1 J g(k)dk (see Section 3.2), which ensures the correct 

o 

behavior of the resulting potential in this region. 

One of our goals was to promote Krein method 4]|S] to solve the inverse problem. 
This method is especially suitable to ascertain the auxiliary potential Vq (r) with no 
bound states, starting from the known Jost function. Having ascertained Vo(r), one 
can build up a series of auxiliary potentials: V\(r), V^(r),..., V n (r), introducing, 
one by one, all bound states with known eigenvalues. Several possibilities of solving 
the main integral equation, Eq. (27), have been proposed and discussed, but not 
yet fully exploited. From the computational-technical point of view the problem 
is much more complicated than it might seem at first sight. Indeed, to ascertain, 
for example, the correct position of the most right-side point in Fig. 5, a system 
of 999 equations has been solved. One can imagine that it is not so easy to extend 
calculations to much larger distances than shown in Fig. 5. Perhaps, to bridge 
over these technical difficulties, one can take some advantage of Eqs. (54) and (56), 
which seem to be straightforward and useful solution schemes. On the other hand, 
rapid development of parallel computing and Grid technology, as well as prospects 
of quantum computing, also suggest some optimism for further research in this field. 
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Figure captions 

Fig. 1. Abraham-Moses [7] trick applied to three-component model potential for 
Xe2 (in ground electronic state) . Each step consists of removing the zeroth 
level along with calculating the new potential according to Eq. (21). The 
depth of the potential well is about 24.3 meV and the original reference 
potential (curve 1) has 24 levels. Although only 5 lowest partner potentials 
(all having the same spectral density for positive energies) are shown, they 
roughly demonstrate that all partner potentials should coincide as r — > 0. 
In addition, the presented curves might help to imagine how the auxiliary 
potential with no bound states would look like. 

Fig. 2. Demonstration of the characteristic function g(k) given by Eq. (5). The 
overall g(k) curve shown in the main figure seems to decay rapidly, but 
this impression is deceptive, because the k scale is logarithmic. As can be 
seen, g(k) = —1 (with high accuracy), if k < k = y/V(0)/C. The insets 
demonstrate how nicely the different intervals can be described by a sum of 
Gaussians (upper inset) or exponents (lower inset) according to Eqs. (30) 
and (35), respectively. Parameters of the different components are given in 
Table 1, and Eq. (39) has been used for the region k > k a = 75000 A. 

Fig. 3. Calculated Krein ii-function in the immediate vicinity of the zero point 
r = 0. 

Fig. 4. Another demonstration of the Krein //-function. The upper graph starts 
where Fig. 3 ends, while the lower graph starts where the upper one ends. 
Note that the period of oscillations is nearly constant and very close to the 

27T 

characteristic value L — — (see Eq. (31)). 

Fig. 5. Calculated G- function as a solution of Eq. (50) (solid curve), i.e., cor- 
responding to the auxiliary potential Vo(r) with no bound states. Note 
that the G-curve falls much slower than the if-curve ascends (cf. with 
Fig. 3). Another curve (open circles) has been calculated according to Eq. 
(59), and is therefore directly related to the reference potential V(r). Since 
these differently calculated curves practically coincide, the two potentials 
in question also coincide in the range depicted. 
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Table 1. Fitting parameters for the characteristic function g(k) 
defined by Eq. (5). In all cases k a denotes the starting point of 
the distance range, while kb marks its end point. 



Mi/A) 


Mi/A) 


3 


cij (dimensionless) 


bj 


Mi/A) 


o 


19230 


1 


0.0642222322635 


28.571730018198 


19266.4518 


based on 


Eq (30) 


2 


0.0427640496056 


31.693733469307 


19214.4537 


bj unit: 


l/A 


3 


1.313714415*10" 3 


18.093176481161 


19169.6355 






4 


-2.18368077*10~ 5 


8.5288122450848 


19153.495 


19230 


20000 


1 


0.0221286191169 


0.0130759575537 


19230 


based on 


Eq. (35) 


2 


0.1065231936499 


0.0024849989851 


19230 


bj unit: 


A 


3 


0.8049658419302 


1.764449543*10" 4 


19230 


20000 


21890 


1 


6.56300549*10- 3 


0.0028620459089 


20000 


based on 


Eq. (35) 


2 


0.1165529570422 


7.050434325*10- 4 


20000 


bj unit: 


A 


3 


0.5953096962391 


1.049812095*10" 4 


20000 


21890 


43600 


1 


0.0384392174818 


5.707750043*10" 4 


21890 


based on 


Eq. (35) 


2 


0.2498847654991 


1.654292141*10" 4 


21890 


bj unit: 


A 


3 


0.2306247842633 


4.074037076*10" 5 


21890 


43600 


75000 


1 


0.0113479356044 


1.341135741*10- 4 


43600 
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Eq. (35) 


2 


0.0540153535949 


5.355103398*10- 5 


43600 


bj unit: 


A 


3 


0.0367554032472 


1.488005999*10" 5 


43600 


75000 


oo 


1 




184142158.24434 
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Eq. (39) 


2 




1695416626*10 16 


(A- 4 ) 






3 




-5.883044*10 24 


(A" 6 ) 
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